An exact solution for 2+1 dimensional critical collapse 
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We find an exact solution in closed form for the critical collapse of a scalar field with cosmological 
constant in 2+1 dimensions. This solution agrees with the numerical simulations done by Pretorius 
and Choputik ||^ of this system. 
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I. INTRODUCTION 

Critical gravitational collapse at the threshold of black 
hole formation, as first found by Choptuik, j2| has been 
studied in many systems. |^] With the exception of a 
study of vacuum, axisymmetric collapse, Q the systems 
studied are spherically symmetric. Because of this sym- 
metry, the equations describing the collapse are PDEs for 
functions of two variables. The critical solution is often 
discretely self-similar (DSS) or continuously self-similar 
(CSS). In the CSS case one can study the critical solution 
itself by assuming the CSS symmetry and thus reducing 
the collapse equations to a set of ODEs. In general, the 
equations (both the PDEs describing collapse and the 
ODEs describing a CSS critical solution) are sufhciently 
complicated that a numerical treatment is needed. 

The Einstein equations in 2-M dimensions are simpler 
than in 3-1-1 dimensions. Though to form a black hole in 
2-1-1 dimensions, one must add a cosmological constant, 
this added feature still allows analytic treatment of some 
collapse situations. [|| One might then hope that critical 
collapse in 2-1-1 dimensions would be more tractable. In- 
deed, for the collapse of thin dust rings or the collision 
of point particles ||^,^ the collapse can be treated ana- 
lytically. (This is essentially because the spacetime has 
constant curvature outside of the zero thickness sources). 

Recently, Pretorius and Choptuik ||l| performed nu- 
merical simulations of the collapse of a massless, mini- 
mally coupled scalar field with a cosmological constant 
in 2-1-1 dimensions. They find that the critical solution 
is CSS. 

In this paper, we find the critical solution of reference 
101 in closed form. In section 2 we write the Einstein- 
scalar equations in an appropriately chosen coordinate 
system. In section 3 we make a CSS ansatz and find the 
solution. This solution is compared in section 4 to the 
numerical results of reference Q and perturbations of the 
solution are considered in section 5. 



II. FIELD EQUATIONS 

The Einstein-scalar equations with cosmological con- 
stant are 



Gab + ^9ab = liTa, 



(1) 



Here, Gab is the Einstein tensor, Tab is the stress energy 
of the scalar field 6 



Ta. 



^a<l)'^b 



1. 



'Jab 



(2) 



and K is a constant. Following the conventions of refer- 
ence [|l| we chose units such that k — 47r. 

We now consider an appropriate choice of coordinate 
system for the metric. Since we want to study the CSS 
critical solution, we want a coordinate system in which 
the solution appears manifestly CSS. This is not the case 
for the coordinates of reference |Q] . Instead we use the 
method of Christodoulou Q to choose a coordinate sys- 
tem where the coordinates are geometric quantities. We 
choose as a radial coordinate f such that 27rf is the length 
of the circles of symmetry. This is the analog of the usual 
area coordinate used in spherical symmetry. We choose 
as our "time" coordinate, the null coordinate u defined 
as follows: u is constant on outgoing light rays, and on 
the world line of the central observer, u is equal to the 
proper time of that observer. Finally, we choose a coor- 
dinate 6 so that d/dO is the Killing vector. The metric 
then takes the form 



(3) 



where v and A are functions of u and r. In reference 
JlOf the Einstein-scalar equations were found for spher- 
ical symmetry in any number of dimensions. For our 
purposes, we specialize the results of reference [0 to 
2-f 1 dimensions, generalize them to add a cosmological 
constant and change the convention to k = 47r. We begin 
by introducing null vectors and ri° defined by 



or 
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du 



dr 



(4) 



(5) 



Then the Einstein-scalar equations with cosmological 
constant are satisfied provided that the scalar field satis- 
fies the wave equation and that the following components 
of Einstein's equation are satisfied 



1 



(6) 



which leads to a scalar field given by 



Gabion'' = A 



Equations and (|^) become 



r or \ or 



-2A " 
2f 9f 



e -■ — (i^- A) = A 



(7) 



(8) 



(9) 



The solution of equations (||) and (^) is most easily ex- 
pressed by defining the quantities g = e"^^ and g = 



. Then we have 



9 = exp 



47r / r\ — ^ I dr 

Or 



9 



1 - 2A / rgdr 
'o 



The wave equation for (f> becomes 
^ d^(t> Ids Id 



dudr r du 



r dr \ ^ dr 



^ 



(10) 



(11) 



(12) 



III. CRITICAL SOLUTION 

We now make the ansatz that the scalar field is CSS 
and use this to solve the field equations. Choose the 
origin of u to be at the singularity, and define two new 
coordinates T and R by 



r = e-'R 
Then demand that (j) take the form 
(j) = cT + V(^) 



(13) 



(14) 



(15) 



where c is a constant. This ansatz requires that we ne- 
glect the cosmological constant, which in turn means that 
g — I and thus reduces equation (12) to the flat space 
wave equation. Putting the ansatz in equation ( |l5|) into 
equation (n2h we obtain 



i?(l - 2R)'4>" -I- (1 - m)^!' 







(16) 



where a prime denotes derivative with respect to R. The 
solution of this equation that is regular at the origin is 



-2c In 



1-F VI -2i? 



(17) 



T - 2 In 



1-1- Vl-2i? 



(18) 



Then using equation ( |10[ ) we find that the metric function 
9 is 



9 = exp 



4W^ R[j^\dR 



(iwr~2^ 

4V1 - 2i? 



(19) 



The metric of the CSS critical solution must be singu- 
lar at u = 0. However, our critical solution (equations 
( p8| ) and (19)) appears to have an additional singularity 
at i? = 1/2 which is the past light cone of the m = singu- 
larity. We now consider whether the apparent singularity 
at i? = 1/2 is a real singularity or a coordinate singular- 
ity. Note that from equation ( p9[ ) it follows that i? = 1/2 
(for any value of T) is a marginally outer trapped sur- 
face and that the Christodoulou coordinates go bad at 
just such surfaces. Now define a new coordinate v by 
V = ~{u -\- 2f). Then the metric is 



H 4irc 



dvdu 



+ ^{u + vfde^ (20) 

Now define the number q and the coordinate w by 
1/(2(7) = 1 ~ 47rc^ and w'^'^ = v. Then the metric is 



ds' = 



—u 



41 l-(2<?)" 



2qdudw 



'de' 



(21) 




This metric is smooth at it; = provided that q — n 
where n is a positive integer. That is, the metric is 
smooth for values of c given by 



(22) 



For n = 1 the spacetime is the 2-1-1 dimensional 
Robertson- Walker metric. 

We now consider the question of whether it is physi- 
cally necessary to impose the condition that the metric 
be smooth at w = 0. If q ^ then one can show using 
equation ( pT]) that the spheres of symmetry are trapped 
surfaces for w < 0. However, the critical solution can- 
not have trapped surfaces, since it forms the boundary 
between those evolutions that result in trapped surfaces 
and those that do not. Therefore, one should expect the 
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numerical critical solution to approach our CSS solution 
only inside the past light cone of the singularity (that is 
for w > 0). It is therefore not physically necessary that 
the CSS solution be smoothly extendible past w — since 
such an extension cannot correspond to the behavior of 
the numerical critical solution. 



IV. COMPARISON WITH NUMERICAL 
RESULTS 

In a near-critical collapse, the evolution at first ap- 
proaches the critical solution, and then diverges from it 
as the single unstable mode grows. Therefore, in compar- 
ing a numerical simulation of near-critical collapse to a 
proposed critical solution, one should make the compar- 
ison at an intermediate time: late enough for the critical 
solution to be approached, but early enough so that the 
unstable mode does not have appreciable amplitude. 

To compare the analytic and numerical results, one 
must express both in the same coordinate system. In the 
coordinates of reference the metric takes the form 



^2A 



{dr^ - d^) + lh&n\r/l)e^^de'^ (23) 



cos2(r/Z) 



where I — —1/ K and A and B are functions of r and 
t. The coordinates of reference |l[ are related to ours by 
f — Ztan(r//)e^ and u is that function oi t ~ r that at 
r = is equal to the proper time of the central observer. 
Therefore, it is fairly straightforward to take scalars and 
tensors in the coordinates of reference |^ and express 
them in our coordinate system. 




FIG. 1. numerical values of d(j)/dT\^^Q (dots) with con- 
stant (line) approximation to intermediate time behavior 



Our solution has an unknown parameter c, which 
should be chosen for best fit with the numerical results. 
To make this choice, we note that the analytical solu- 
tion (equation (p^) satisfies c — dcj)/ dT\^^f^. Thus the 
values of 90/9T|^^q for the numerical solution should 
allow us to determine the value of c. Figure 1 shows a 
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FIG. 2. comparison of analytic (curve) and numerical 
(dots) values of the scalar field at an intermediate time 



plot of d(j)/dT\j^^^ vs T for the numerical solution. Note 
that there is a range of intermediate times for which this 
quantity is approximately constant. The line in the solu- 
tion corresponds to the value of c given by equation ( [2^ ) 
with n ~ 4. While it is not necessary that c be given 
by equation (p^), we find that this particular value of c 
gives excellent agreement with the numerical data. From 
now on, we will assume that c has this value. 

Figure 2 shows a comparison between the numerical 
and analytic results for the scalar field (p. Here the dots 
are the numerical results and the curve is the analytic 
one. The comparison is made at T = 9. The freedom 
to add a constant to the scalar field is used to set the 
value of the scalar field to zero at the origin. Note that 
there is excellent agreement between analytic and numer- 
ical results. Though this figure shows the comparison at 
only one time, the agreement persists for a large range of 
intermediate times. 

Due to the symmetry of the spactime, the metric is de- 
termined entirely by the matter. This is made explicit in 
equations ( p^ and (ll|). Therefore, since the scalar fields 
of the numerical and analytic solutions agree, the metrics 
must agree also. Nonetheless, for illustrative purposes we 
present a comparison of the metrics of the numerical and 
analytic solutions. Figure 3 shows a comparison of the 
quantity M = {f'^/P) — [g/g) with the numerical solu- 
tion given by the dots and the analytic solution given by 
the curve. Note the excellent agreement between the two 
solutions. Here the comparison is made at T = 9, but 
the agreement persists for a large range of intermediate 
times. The metric is determined by M and g which is 
1 for the analytic solution and very close to 1 for the 
numerical solution in the intermediate range of times. 
Thus, in the intermediate range of times, the scalar field 
and metric of the numerical and analytic solutions agree. 

In figures 2 and 3, we have performed the compar- 
ison of analytic and numerical solutions only up to the 



3 




FIG. 3. comparison of analytic (curve) and numerical 
(dots) values of the mass aspect at an intermediate time 



past light cone of the singularity. However, the numerical 
solution certainly continues beyond this light cone, and 
since q = 4, the analytic solution also extends. Do the 
two solutions still agree in this region? I argued at the 
end of the previous section that they cannot agree, since 
the analytic solution contains a closed trapped surface 
beyond the light cone. Nonetheless, it would be useful to 
compare the two solutions in this larger region to see the 
extent and nature of their disagreement. 

To make this comparison, a different coordinate from 
R is needed, since the solution is not smooth in R across 
the light cone. Here we use as our new coordinate the 
affine parameter A along an outgoing radial null geodesic. 
The affine parameter is made scale invariant by choosing 
the null geodesic to have inner product u with the central 
observer. The analytic solution can be expressed para- 
metrically in terms of A as follows: (here we specialize to 
the case q = 4). Introduce the variable x = {—u) ' w. 
Then x — 1 at the origin and a: = at the past light cone 
of the singularity. We have 



A 



1 



2V2 



1 



x^V^^dx 




M 



1 



7/2 



(24) 



(25) 



(26) 



From equation ( p4| ) it follows that the past light cone is 
at A w 0.838. Using equation (p3|), the affine parame- 
ter A can be found numerically from the numerical data. 
Note, however that the numerical simulations have only 
a certain range in A since a critical numerical solution 



must stop before the time t of singularity formation. In 
our case, the maximum value of A for the numerical data 
is approximately 1. 

Figures 4 and 5 show a comparison between the nu- 
merical and analytic solutions for the range < A < 1. 
Here the values of 4> are compared in figure 4, while the 
mass aspect M is compared in figure 5. As in figures 2 
and 3, the comparison is made at T = 9. From the figures 
it is not clear whether there is a disagreement between 
the two solutions in the range A > 0.838 that is beyond 
the past light cone. To see the disagreement clearly, one 
would need a larger range of A. This could perhaps be 
provided by a numerical simulation that used double null 
coordinates. 




FIG. 4. comparison of analytic (curve) and numerical 
(dots) values of the scalar field with respect to affine param- 
eter at an intermediate time 




FIG. 5. comparison of analytic (curve) and numerical 
(dots) values of the mass aspect with respect to affine pa- 
rameter at an intermediate time 
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V. PERTURBATIONS 

We now turn to a treatment of perturbations of the 
critical solution. This treatment should in principle allow 
us to determine analytically (1) the correct value of the 
parameter c and (2) the value of the exponent in scaling 
relations for near-critical collapse. Unfortunately, as we 
will see below, it is unclear what boundary conditions to 
impose on the perturbations. 

The critical solution, when perturbed, has one unstable 
mode that grows as e'^'-^ for some constant k. Therefore 
the correct value of c is the one for which there is exactly 
one unstable mode. The quantity k is related to scaling 
laws for near-critical collapse. In the collapse of a one pa- 
rameter family of initial data, a quantity Q with dimen- 
sion (length)'* obeys a scahng relation Q cx \p — p*\^^'' 
where p is the parameter and p* is its critical value. 
Pretorius and Choptuik examine scaling in maximum 
scalar curvature for subcritical collapse, a quantity that 
has dimension (length)"^. They find k « 0.83. 

In the approximation that g = 1, the scalar field sat- 
isfies the flat space wave equation. Therefore, the per- 
turbed scalar field S(p also satisfies this equation. Then 
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making the ansatz ( 
we obtain 



e S{R) and using equation (113) 



R{1 - 2R)S" + {l-[3 + 2k]R)S' - kS ^ (27) 
The solution is 



5 = i^(fc,l/2,l,2i?) (28) 

where -F is a hypergeometric function. Writing the hy- 
pergeometric function in integral form, we have 

5 = - / (1 - 2Rsin^x) dx (29) 

71" Jo 

In principle, one should now impose boundary condi- 
tions on the behavior of the perturbation at i? = 1/2 and 
these boundary conditions would then determine the al- 
lowed values of k. Unfortunately, it is not clear what 
boundary conditions are physically reasonable, since the 
critical numerical solution should match the analytical 
solution only for R < 1/2. Therefore, though we have 
the solution of the perturbation equation, we cannot use 
this solution to determine the critical exponent. 
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